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A  CLASSICAL  ELECTROMAGNETIC  APPROACH 


by  Daniel  A.  Jelski,  P.  T.  Leung  and  Thomas  F.  George 


Departments  of  Physics  &  Astronomy  and  Chemistry 
239  Fronczak  Hall 

State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260  USA 


ABSTRACT 

This  review  article  discusses  several  important  aspects  of  photochemistry  a 
structured  metallic  surfaces.  The  electromagnetic  field  above  the  surface 
is  calculated  using  the  Rayleigh  expansion.  Conditions  under  which  this 
expansion  is  valid  and  simplifying  approximations  which  make  it  easier  to 
use  are  discussed  in  detail.  This  formalism  is  then  applied  to  three 
different  but  related  phenomena.  First,  the  photodissociation  rate  of  a 
molecule  above  a  surface  is  calculated  for  laser  frequencies  at  or  near  the 
surface  plasmon  resonance.  It  is  found  that  there  is  an  optimal  molecule- 
surface  distance  for  photodissociation.  Then  the  absorption  lineshape  of  a 
molecule  is  considered,  where  both  Fano  and  Lorentzian  lineshapes  are  found 
to  be  distorted  as  the  molecule  approaches  the  surface.  Finally,  laser- 
induced  periodic  deposition  is  discussed,  and  a  model  is  developed  to 
describe  the  growth  rate  of  a  cadmium  grating. 


1. 


Introduction 


Catalysis  and  control  of  chemical  reactions  is  one  of  the  most 
important  tasks  of  a  chemist.  Many  different  methods  have  been  tried  to 
achieve  both  selectivity  and  high  product  yields.  One  of  the  best  methods 
now  available  is  laser  chemistry,  which  permits  the  chemist  a  great  deal  of 
control.  Laser-controlled  chemical  reactions  have  been  studied  for  the  past 
decade  (George,  1982).  A  second  popular  method  of  catalyzing  reactions  is 
on  surfaces,  in  particular  metal  surfaces  such  as  platinum  have  long  been 
used  in  organic  chemistry.  It  is  thus  timely  and  important  to  inquire  into 
the  dual  effects  of  these  mechanisms  combined.  Hence  the  topic  of  this 
review  article  concerns  the  interaction  of  laser  light  with  a  surface,  and 
the  chemistry  which  can  thereby  be  accomplished. 

While  more  conventional  laser  chemistry  has  long  been  done  in  bulk  and 
gas  phase  material,  the  possibility  of  enhancing  reaction  rates  near  a 
surface  has  been  a  tantalizing  prospect  ever  since  the  discovery  of  the 
dramatic  surface-enhanced  Raman  scattering  (SERS)  (Fleischmann  et  al.,  1974) 
from  pyridine  adsorbed  on  silver  electrodes.  It  has  subsequently  been 
learned  that  the  most  important  single  factor,  at  least  for  physisorbed 
molecules,  in  achieving  this  surface  enhancement  is  the  roughness  of  the 
substrate  surface.  Surface  roughness  couples  with  the  radiation  field  to 
produce  a  plasmon  along  the  surface,  which  in  turn  produces  a  large  field 
near  the  surface  under  resonance  conditions.  This  accounts  for  the  SERS 
effect  just  described. 

The  possibility  of  using  this  enhanced  surface  field  to  effect 
photochemistry  above  and  on  the  surface  has  been  much  studied  for  ten  years. 
A  review  of  the  experimental  literature  was  written  by  Goncher  et  al. 


(1984).  The  present  review  concentrates  on  some  theoretical  aspects  of  the 
problem,  and  discusses  some  models  to  account  for  photodissociation  above  a 
surface  and  laser-induced  deposition  on  a  surface. 


The  major  concern  here  is  the  application  of  classical  electromagnetic 
(EM)  theory  to  the  study  of  possible  photochemical  processes  occuring  in  the 
vicinity  of  a  structured  metallic  surface.  To  this  end,  we  and  others  have 
found  the  Rayleigh  expansion  of  the  surface  EM  fields  invaluable  in  making 
the  calculations  both  tractable  and  accurate.  In  the  next  section  we  review 
when  and  under  what  circustances  the  Rayleigh  expansion  is  valid  and  shall 
also  discuss  various  approximations  which  make  practical  calculations 
easier.  In  Sec.  3  we  consider  an  application  of  the  expansion  to  the 
calculation  of  the  photodissociation  dynamics  of  a  molecule  above  a  rough 
surface.  In  this  case  our  use  of  the  Rayleigh  expansion  is  justified 
everywhere  as  we  are  considering  a  shallow  grating.  In  Sec.  4  we  consider 
the  absorption  lineshape  of  a  molecule  above  a  surface.  Here  we  also 
consider  the  coupling  between  the  molecule  and  the  surface  plasmon  along  a 
shallow  grating.  In  Sec.  5  we  use  the  Rayleigh  expansion  where  it  is  exact 
(above  the  selvedge  region)  but  also  employ  an  approximate  method  of 
calculation  valid  for  deeper  gratings.  This  will  permit  us  to  discuss 
laser- induced  chemical  vapor  deposition  as  recently  observed  (Brueck  et  al., 
1982).  Finally,  there  is  a  brief  conclusion  in  Sec.  6. 
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2 .  The  Rayleigh  Expansion 

The  Rayleigh  expansion  has  been  widely  used  and  discussed  in  the 
literature  (Rayleigh,  1907;  Fano,  1941;  Petit,  1980).  Briefly,  this 
expansion  is  a  solution  of  the  homogenous  Helmholtz  equation 


Au  +  k  u  =  0  , 


(2.1) 


where  u  represents  the  electric  field  £  or  the  magnetic  field  S,  depending 
on  the  polarization.  Maxwell's  equations  reduce  to  this  simple  form  for 
time -harmonic  fields  in  the  absence  of  sources.  Hence  for  light  irradiating 
a  rough  surface,  the  Rayleigh  expansion  is  an  exact  solution  to  Maxwell's 
equations  outside  the  selvedge  region.  Inside  the  selvedge  region  the 
situation  is  much  more  complicated  given  the  charges/currents  along  the 
surface.  Then  Eq.  (2.1)  becomes  inhomogenous  and  other  methods  have  to  be 
found  to  solve  the  equation. 

The  Rayleigh  hypothesis  states  that  the  Rayleigh  expansion,  exact  above 
the  selvedge  region,  is  also  a  solution  within  the  selvedge  region  as  long 
as  the  grating  is  shallow.  This  turns  out  to  be  reasonably  accurate  (Glass 
et  al.,  1981)  provided  E,^/ A  <  0.1,  where  E ^  is  the  grating  amplitude  (more 
generally,  we  define  En  as  being  the  Fourier  amplitude  of  the  k^  wavenumber) 
and  A  is  the  grating  wavelength.  Beyond  this  point  the  Fresnel  matrix 
becomes  ill-conditioned  and  hence  precludes  a  solution  of  the  Rayleigh 
equations.  Agassi  and  George  ( 1986a, b)  have  developed  a  means  to  prevent 
tliese  numerical  difficulties  and  are  able  to  solve  the  Rayleigh  equations 
for  gratings  of  arbitrary  height.  They  go  on  to  maintain  that  this  is  in 
fact  the  exact  solution  for  certain  kinds  of  gratings,  even  within  the 


selvedge  region.  While  we  now  begin  to  question  this  latter  statement,  we 
shall  nevertheless  find  Agassi  and  George's  dressed  expansion  to  be  very 
useful  outside  the  selvedge  region  where  its  validity  is  guaranteed. 

Let  u  of  Eq.  (2.1)  refer  to  the  electric  field  (p-polarization) . 
Suppose  also  that  we  have  a  periodic  grating  with  wavenumber  k^.  The 
incident  light  with  wavenumber  k  =  Wq/c  is  of  intensity  E^.  The  surface 
lies  along  the  z  =  0  plane,  separating  two  dielectrics  of  complex  constants 
e(0)  and  e(l)  above  and  below  the  surface  respectively.  Without  loss  of 
generality  we  shall  henceforth  assume  that  e(0)  =  1,  and  hence  we  can  refer 
to  e(l)  simply  as  e.  Then  we  can  write 


-a  =  c0(0)?a.(0)exp(i[kox-V))  +  ^  A£pa+U)exp(  i[kjix+aiz] ) 

(2.2) 


l=-° 


z  >  £(x) 


Ib  =  \  Pg+UKjexpUlkjX-S^J) 


Z  <  s(x), 


where 


k  =  ^ 


el  +  iE2 


<  -1 


ke  -  kt  +  %  \  *  (fc2  ’  kt2)1/2 


Ima  >  0 
a„  >  0 


(2.3) 
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2  2  1/2 

^  =  (ek  -  )  ImS^  >  0 

pa±u)  =  (ke£  ;  V1  P3±(°  =  k[k^  ;  M3  • 

A  similar  pair  of  equations  exists  for  the  magnetic  field. 

Aq(2.)  and  C^(i)  can  be  found  by  matching  the  boundary  conditions.  It 
is  clear  that  Cg(0)  represents  the  amplitude  of  incoming  waves  and  hence 

Cn(0)  =  E.  .  (2.4) 

0  l 


Toigo  et  al.  (1977)  have  devised  an  elegant  method  for  solving  the  coupled 
equations  that  result  when  Eq.  (2.2)  is  matched  across  the  boundary.  The 
derivation  is  repeated  in  Agassi,  et  al.  (1986a).  This  solution  requires 
evaluating  integrals  of  the  following  form 


1  fX,  i  ( a 
t  dx  e  n 

A  J0 


-3  )5(x)  -i(k  -k  )x 


m 


n  rn  i  $(n-m) 


(2.5) 


where  £(x)  is  the  surface  profile  function.  This  integral  is  soluble 
analytically  for  £(x)  a  sinusoidal  or  sawtooth  function  (Laks,  1981). 

The  solution  for  the  sinusoidal  case,  where  is  the  amplitude  of  the 
grating,  is  given  by  Agassi  and  George  as 


V  M  .A.  =  y  E . , 
m,£  l  mi 

t«-oo 

(2.6) 


%  %  » 


y  « 

L  m> 
e=-« 


lkl  =  VmEi’ 


where 


aoe  +  k.k 

^  ®  2-  01  /  ,  \JD*i  «  /  -  r  ^  i  \ 

Mm  9  ~o -  (i)  j  0U,  [a  -  8  ]) 

m,  t  a  -  p  m-t  lx.  m 


m 


a.  8  +  k.k 

ii  2.  in  ?.  m  /  .  \in“  9. ,  / ,.  r  \ 

N  9  - -  U)  J  *  M) 

p_ot  iq*  t,  1  m  £ 


m 


-ang  +  k-.k 

u  m  um/1-vrn*/»-r  »-*  i\ 

“m  ‘  an  4-  |S  <l)  J(5llt,0  +  BJ> 
0  m 


(2.7) 


v 

m 


2c 

e-1 


a  6  n 
m  m,0 


The  coefficients  A^  and  can  be  determined  by  inverting  the  matrix 
equations  (2.6). 

Agassi  and  George  have  developed  a  scheme  to  assure  convergence  of  the 
Rayleigh  expansion  under  all  circumstances.  Toigc  et  al.  (1977)  have  shown 
that  it  should  be  possible  to  do  this.  Agassi  and  George  (1986)  point  out 
that 


Ateiai5(x)  ^  Aie'kglilC(x)  for|l|->-  .  (2.8) 

Now  suppose  that  £(x)  has  a  large  amplitude.  Then  when  £(x)  <<  0  the 
exponential  diverges,  and  to  ensure  the  convergence  of  the  total  field,  A^ 


must  get  exponentially  small.  This  implies  that,  in  order  to  calculate  the 
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Rayleigh  coefficients,  we  must  invert  a  matrix  where  elements  M  are 

m,  x. 

growing  exponentially  in  size.  Clearly  this  becomes  numerically  unstable. 
We  can  correct  this  deficiency  by  making  the  transformation 


=  A^e"ia2.^£  and  =  C^e"16^^  . 


We  can  then  transform  Eqs.  (2.7)  as 


(2.9) 


M  =  M.  e^am^m  and  N  =  N  e^m^m  , 
£,m  £,m  i,m  i,m 


(2.10) 


whereby  Eq.  (2.6)  becomes 


y  m 

L  m* 


.A  =  u  E. 

I  l  mi 


(2.11) 


which  converges  for  arbitrary  grating  depth. 

An  alternative  form  for  the  Rayleigh  expansion  can  be  derived  using 
Green's  functions.  The  advantage  here  is  that  one  also  obtains  a  general 
form  for  the  solution  to  the  inhomogenous  equation  which  reduces  to  the 


Rayleigh  expansion  outside  the  selvedge  region,  as  demonstrated  by  Toigo  et 
al.  (1977).  The  form  for  this  equation  outside  the  selvedge  region  is 


m  k“  -  k.k 

- a — ”  H(n)  "  iL(m)]*U-»)  =  2alEi6t,( 


(2.12a) 


v  *.•  v  v  v  \* 


%  s 


H(m)  -  iL(m)  ] 4>(  H-m)  =  0 


(2.12b) 


The  disadvantage  of  this  formalism  is  that  it  involves  two  coupled  equations 
which  are  not  separable  in  the  way  that  equations  (2.6)  are.  One  must  solve 
for  both  H(m)  and  L(m)  simultaneously. 

Recently  a  perturbative  technique  has  been  developed  using  Eq.  (2.12) 
which  permits  the  solution  of  field  intensities  in  the  shallow  grating 
limit  (Weber,  1986).  The  primary  purpose  was  to  calculate  the  dispersion 
relation  for  surface  plasmons  along  a  rough  surface.  Weber  (1986)  has 
extended  this  idea  to  an  approximate  (but  quite  accurate)  calculation  of  the 
reflectivity.  From  his  paper  we  find  we  can  write 


m)  =  y  - n  m  [A1 . (m,n)L(n)  +  A._(m,n)H(n)J  - 


ea  +p 
n  mm 


2ia  k  H( 0) 
u  z _ 

eao+po 


6m  n  (2-13) 
m,  U 


r  l  2ia  k  H(0) 

=  \ - !1ISL.[A  (mfI1)L(n)  +  A__(m,n)H(n)  ] - — -  6  n>  (2.14) 

L  +3  21  22  ean+3n  m’° 

n  m  m  0*0 


where  |  is  a  function  of  <$(m,n)  which  reduces  to  the  (n-m)-th  Fourier 
n-m 

component  of  £(x)  in  the  shallow  grating  limit.  A,,  is  a  2  *  2  matrix  which 
depends  only  on  the  quantities  defined  in  Eq.  (2.3),  and  kz  is  the 
perpendicular  component  of  the  incident  field.  As  mentioned  previously, 


these  equations  are  fundamentally  the  same  as  equations  (2.6)  above  the 
selvedge  region. 

Let  us  now  define  the  vector 


V  5  >!?»i 

m 


(2.15) 


whereby  we  can  write  Eqs.  (1.14)  and  (1.15)  as 


V  =  <t>0 6  .  +  )  ,  _  f 

m  0  m,0  /  2  2  n 

u  w  -  OJ 

n  m 


(2.16) 


where 


.  2.  2  e+1 

to  =  c  k  - 

m  me 


and  is  a  vector  containing  the  rightmost  terms  of  Eqs.  (2.14)  and  (2.15). 

V  is  a  2  *  2  matrix  which  depends  on  A  and  also  on  £ 
mn  r.-m 

Now  suppose  that  three  modes  are  at  or  near  resonance  with  the  surface 
plasmon.  Obviously  we  have  to  include  the  specular  reflection  as  being 
significant  since  we  are  interested  in  calculating  the  reflectivity. 
Furthermore,  we  should  include  two  resonant  frequencies,  for  at  least  +kr 
and  -k  will  be  resonant.  We  can  then  separate  out  these  terms  from  Eq. 


(2.16)  to  yield 


4,0<5m,0  + 


n=0 , ra, , m.  m 


+  )  -0-  mn  -  V 

/  2  2  n 

“  u>  -  U) 

n*0,m, ,m_  m 


(2.17) 


For  non-resonant  terms  (n?=0(m^,mo)  we  can  make  the  approximation 


¥ 

P 


l 

n-0, m  ^ f 


V 

■Pn  - 

2  2 
oj  -  UJ 

P 


p  *  0,  m^ ,  m0 


(2.18) 


Substituting  this  back  into  Eq.  (2.17)  and  doing  some  algebra,  we  get 


^m  *06m,0  + 


[V  + 
mn 


n “0  f  m  ^  >  m0 


■■^P  P~  ]¥  (ld2  -  -d2)'1/2.  (2.19) 
2  2  n  m 

_  OJ  -  U1 

p»0,m1,m2  p 


Recalling  that  is  a  2  x  2  matrix,  we  are  led  to  a  6  *  6  matrix  equation 

for  ¥  .  However,  for  modes  that  are  resonant  we  can  use  the  relation 
m 

3 

L(n)  *  —  H(n)  for  n  =  0,  m^,  •  (2.20) 

Then  the  system  reduces  to  a  3  *  3  matrix  which  can  be  solved  analytically. 

Weber  numerically  tested  this  method  against  the  exact  theory.  His 
model  was  a  silver  surface  with  a  sawtooth  grating  profile,  A  =  8000  A.  Th 
incident  frequency  was  varied  between  about  2.15  and  2.25  eV,  with  the 
"plasmon  dip"  in  the  reflectivity  observed  at  about  2.21  eV  (see  Fig.  1). 
For  a  grating  height  of  300  A  the  approximate  method  corresponds  almost 
exactly  to  the  exact  calculation.  For  =  300  A  some  error  is  introduced 
but  the  method  still  works  remarkably  well.  The  result  as  a  function  of 


grating  height  is  shown  in  Table  1. 


Our  group  has  found  an  even  simpler  method  for  calculating  the  plasmon 

field  intensity  at  resonance  (Jelski  et  al.t  1987).  Rather  than  use  the 

extinction  theorem  results  of  Eq.  (2.12)  as  Weber  did,  we  instead  use  the 

Rayleigh-Fano  expansion  of  Eq.  (2.2).  This  implies  that  we  need  not 

introduce  the  vector  as  Weber  does,  but  since  the  equations  are  separable 

n 

we  obtain  a  scalar  equation.  However,  since  we  restrict  ourselves  to  the 
resonance  frequency,  we  can  further  simplify  the  problem  by  neglecting 
specular  reflection.  Hence  our  method  will  not  yield  accurate  values  for 
the  reflectivity  away  from  resonance.  We  have  developed  it  to  calculate  the 
plasmon  intensity  as  a  function  of  grating  height.  In  this  application  it 
appears  quite  suitable  (see  Fig.  2),  and  the  results  compare  favorably  with 
those  of  Weber.  We  shall  investigate  the  application  more  in  Sec.  5. 

Unlike  Weber,  our  formalism  is  applied  to  a  sinusoidal  surface,  and  hence  we 
can  use  the  results  embodied  in  Eqs.  (2.6)  and  (2.7). 

We  begin  with  Eqs.  (2.6).  Let  us  assume  for  specificity  that  is 
resonant,  and  hence 

A.  >>  A  m  *  1 .  (2.21 ) 

1  m 

Then  we  can  rewrite  Eqs.  (2.6)  as 


l 

+  L  ,  M, 


M,  E.  +  u  ,  M,  A 
1  1  m*l  lm  m 


(2.23) 


The  similarity  in  spirit  to  Weber's  derivation  should  now  be  obvious. 
Continuing  in  that  vein,  we  approximate  A^  as 


U  E.  -  M  , A. 
mi  m ,  1  1 


and  substituting  this  into  Eq.  (2.23),  we  obtain 


(2.24) 
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(2.24) 


For  normal  incidence,  however,  it  is  not  possible  to  ignore  the  A  ^  term 
since,  by  symmetry,  A^  =  A Including  this  effect  we  get 


<>*i  '  [  <MiP  +  "i-p'H®  >Ei 
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p*l  pp 


(2.26) 


To  ensure  convergence  we  can  use  the  dressed  form  of  these  equations  as 
previously  discussed. 

The  results  of  this  method  are  shown  in  Fig.  2.  We  see  that  the 


approximate  method  overstates  the  exact  result,  but  that  it  is 


qualitatievely  correct.  On  comparison  with  Weber's  data  (Table  I),  we  find 
the  same  qualitative  behavior  despite  the  differences  in  the  model  studied. 
How  does  our  method  compare  with  that  of  Weber?  The  first  obvious  fact  is 
that  we  have  neglected  the  specular  reflection  from  our  calculation.  Hence 
our  calculation  of  the  reflectivity  is  obviously  not  possible.  However,  at 
plasraon  resonance  the  reflection  may  be  considered  small  and  hence 
neglected.  Thus  our  method  is  thus  suitable  for  calculating  resonance 
plasmon  field  intensities.  Weber  has  also  included  two  other  resonant 
frequencies,  and  n^.  We  have  done  the  same  for  the  case  of  normal 
incidence.  For  other  cases  it  would  seem  unlikely  that  there  would  be  two 
resonant  frequencies  very  close  together.  Hence  the  restriction  of  normal 
incidence  is  probably  not  severe  from  the  standpoint  of  prac  :cal 
application.  In  conclusion,  we  note  that  while  Weber's  method  is  more 
general,  ours  is  easier  to  use  and  it  seems  that  the  derivation  is  more 
transparent.  Further,  under  conditions  where  our  method  applies  it  appears 
to  be  just  as  accurate. 

We  continue  our  discussion  of  the  Rayleigh  expansion  by  considering  an 
important  special  case,  namely  the  shallow  grating  limit.  In  this  case 
♦(n-m)  can  be  expanded  to  first  order  in  5(x) .  This  has  been  discussed  in 
an  excellent  review  by  Maradudin  (1982),  and  it  is  also  the  approach  taken 
by  Jha  et  al.  (1980).  We  shall  use  these  results  in  Sec.  3.  To  first  order 
in  £(x),  we  have 


♦(n-m)  =  <5  + 

n-m 
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ra  n-m 


(2.27) 


In  this  case,  Eqs.  (2.7)  simplify  to 
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M  =  im(-a_g  +  knk  )E, 
m  CTm  Omm 


If  we  insert  these  results  into  Eq.  (2.26),  we  get  the  first-order 
approximation  for  the  plasmon  field  intensity.  The  result  (with  different 
notation)  is  given  in  the  next  section  by  Eqs.  (3.6)  and  (3.7)  and  will  be 
discussed  further  in  Sec.  5.  Similarly,  by  setting 


(2.29) 


one  recovers  the  flat  surface  dispersion  relation.  For  real  e  this  is  given 
by  Eq.  (3.8).  The  case  of  complex  z  will  be  discussed  in  Sec.  5. 

We  close  by  considering  some  attempts  to  solve  the  inhomogenous 
equation  valid  in  the  selvedge  region.  This  has  been  done  analytically  for 
the  case  of  the  square  well  grating  (Sheng  et  al.,  1982).  The  result  has 
been  extended  (Lee  et  al.,  1985)  by  taking  advantage  of  the  scheme  depicted 
in  Fig.  3.  It  is  possible  to  determine  the  field  strength  at  each  layer 
(exactly)  from  the  result  for  the  previous  layer.  A  recursion  relation  thus 
develops  which  gives  the  total  field  everywhere  exactly. 

3.  Photodissociation 

While  both  vibrational  and  electronic  molecular  spectroscopy  on  rough 
metallic  surfaces  have  been  studied  extensively  (Avouris  et  al.,  1984; 
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Moskovits,  1985),  the  general  area  of  photochemistry  at  such  surfaces  has 
also  attracted  much  activity  in  recent  years  (Goncher  et  al.,  1984).  In 
particular,  the  process  of  photodissociation  of  gas  molecules  at  such 
surfaces  is  of  great  interest  since  it  is  the  first  step  that  one  must  study 
in  order  to  understand  and  control  the  various  phenomena  ranging  from  the 
deposition  of  molecules  (Brueck  et  al.,  1982)  (see  Sec.  4)  to  laser-induced 
heterogeneous  catalysis  (Lin  et  al.,  1984).  In  this  section,  we  shall 
explore  briefly  both  experimental  and  theoretical  work  done  in  recent  years 
in  this  area.  We  shall  emphasize  the  roughness  and  surface  plasraon 
excitations  of  metallic  surfaces,  in  constrast  to  some  previous  reviews  on 
spectroscopy  where  other  excitations  such  as  electron-hole  pairs  and  phonons 
have  been  stressed,  in  which  cases  the  role  of  the  surface  roughness  has 
been  minor  (Avouris  et  al. ,  1984).  Furthermore,  we  shall  limit  ourselves 
mainly  to  physisorbed  molecules  where  bonding  between  the  admolecule  and  the 
surface  can  be  neglected.  For  the  case  of  chemisorbed  molecules  due  to 
charge-transfer  processes  occurring  between  the  admolecule  and  the  surface 
(Lundqvist,  1984),  the  problem  has  become  more  complicated,  since  under 
these  circumstances  photodesorption  will  accompany  photodissociation  of  the 
admolecules  as  a  competing  process,  which  requires  a  more  difficult 
theoretical  analysis.  Nevertheless,  experimental  work  has  been  done  along 
these  lines  (Bourdon  et  al.,  1984,  1986). 

We  begin  by  reviewing  some  experiments  recently  done  on  the 
photofragmentation  of  physisorbed  molecules  on  rough  metallic  surfaces.  The 
first  one  we  would  like  to  mention  is  that  carried  out  at  Lincoln  Laboratory 
(Ehrlich  et  al.,  1981)  in  which  metal-alkyl  compounds  are  deposited  on  a 
host  substrate  with  both  gas  and  liquid  adlayers  coexisting  on  the 


substrate.  Photodissociation  of  such  absorbates  is  studied  by  using  a  UV 

laser  of  weak  intensity  (-3  mW)  at  257.2  nm  and  for  two  different  metal 

alkyls,  dimethyl  cadmium,  Cd(CH_.)-,  and  hexamethvl  alminium,  A1-,(CH-,),.  It 

J  l  '  l  j  o 

is  found  that  for  the  case  of  CdCCH^)^  ,  most  of  the  dissociation  occurs  in 

the  gas  phase  well  above  the  substrate,  and  for  the  case  of  Al^CCH^)^  most 

of  the  dissociation  occurs  well  into  the  adlayers.  This  seems  to  imply  an 

optimal  molecule-surface  distance  at  which  maximum  dissociation  occurs. 

Indeed,  in  performing  an  experiment  of  photodecomposition  of  pyridine 

molecules  adsorbed  on  roughened  silver  Ag(110)  surface  at  406.7  nm  laser 

frequency,  a  group  at  Berkeley  (Goncher  et  al.,  1984)  has  observed  the 

existence  of  such  an  optimal  distance.  The  next  experiment  we  want  to 

mention  is  that  carried  out  by  the  Exxon  group  (Garoff  et  al.,  1982),  where 

photodegradation  of  dye  molecules  ( rhodamine  6G)  by  a  visible  laser  of  -0.1 

2 

W/cm  has  been  carried  out  on  top  of  a  silver- island  film  on  a  silica 
substrate.  Enhanced  photodegradation  is  hardly  observed  in  this  experiment, 
but  rather  for  molecules  close  to  the  silver  islands,  decreased 
fragmentation  rate  is  seen.  On  the  other  hand,  an  experiment  conducted  by 
Columbia  researchers  has  reported  the  observation  of  enhanced 
photodissociation  of  organometallic  molecules  at  metallic  island  surfaces 
(Chen  et  al.,  1983).  In  this  experiment,  the  same  UV  (257.2  nm)  laser 
(Ehrlich  et  al.,  1981)  is  used  to  dissociate  Cd^H^^  on  top  of  a  dielectric 
substrate  which  is  covered  by  a  mixture  of  spherical  cadmium  and  gold 
pellets.  These  spheres  are  observed  to  grow  to  ellipsoids  due  to  subsequent 
deposition  of  the  Cd  molecules  following  the  dissociation  process.  Enhanced 
growth  rate  has  been  observed  for  Cd  spheres  but  not  for  Au  spheres.  As  we 
shall  see  below,  these  seemingly  contradictory  observations  (Garoff  et  al., 


1982;  Chen  et  al.,  1983)  can  be  explained  by  introducing  the  concept  of  a 
critical  molecule-surface  distance  into  the  description  of  the  dissociation 
phenomena.  Let  us  first  review  briefly  some  theoretical  work  which  helps 
serve  as  a  basis  for  such  a  concept. 

Ever  since  the  first  observation  of  surface-enhanced  Raman  scattering 
(SERS)  (Fleischman  et  al.,  1974),  a  large  amount  of  effort  has  been  devoted 
to  the  theoretical  explanation  of  this  phenomenon  and  the  investigation  of 
the  possible  surface  enhancement  of  other  photochemical  processes  (Chang  et 
al.,  1982;  Moskovits,  1985).  To  this  latter  effort,  we  want  to  mention  in 
particular  the  work  by  Nitzan  et  al.  (1981),  Weitz  et  al.  (1983)  and  Gersten 
et  al.  (1985),  who  have  studied  both  resonant  and  non-resonant  processes 
including  Raman,  resonance  Raman,  fluorescence  and  photoabsorption 
phenomena.  In  most  of  these  cases,  the  four-level  model  has  been  found  to 
be  quite  successful  in  the  explanation  of  these  various  phenomena  (Weitz  et 
al.,  1983).  The  main  physical  mechanisms  have  been  identified  to  include 
the  image  effect,  the  shape  (lightening  rod)  effect  and  the  surface  plasmon 
(SP)  effect.  While  the  first  one  is  found  to  be  very  small  in  ordinary 
SERS,  the  last  one  is  viewed  as  the  main  mechanism  leading  to  such  dramatic 
enhancement.  This  is  true  at  least  for  physisorbed  molecules  in  which  the 
bonding  effect  between  the  admolecule  and  surface  can  be  neglected,  though 
the  situation  may  be  different  for  chemisorbed  molecules  (McCall  et  al., 
1980).  Two  conditions  must  be  met  in  order  for  the  plasmon  enhancement 
effect  to  be  plausible  ,  namely,  the  metallic  surface  must  be  rough  and  the 
incident  light  frequency  should  satisfy  the  SP  resonance  condition. 
Roughness  implies  the  excitability  of  SP,  and  the  resonance  condition 
ensures  a  large  magnitude  of  the  SP  field  which  then  leads  to  strong 


absorption  by  the  molecular  system  on  top  of  the  surface.  These 


investigators  have  modelled  the  surface  roughness  as  a  collection  of  very 

tiny  spneres  and  have  provided  specific  details  for  an  isolated  sphere  and 

the  case  with  two  neighboring  spheres  (Gersten  et  al.,  1985).  Furthermore, 

they  have  extended  their  investigation  to  the  possibility  of  surface- 

enhanced  photodissociation  for  the  case  of  direct  (fast)  dissociation.  This 

latter  process  can  be  treated  on  the  same  footing  as  the  absorption  process 
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since  the  dissociation  takes  place  on  a  time  scale  of  the  order  of  10  s 
following  absorption,  leading  to  a  yield  of  almost  unity  for  such  reactive 
processes  (Nitzan  et  al.,  1981).  To  calculate  this  absorpt ion/dissociation 
cross  section,  the  classical  approach  has  been  taken  in  which  the  molecule 
is  modelled  as  a  point  dipole  (u)  satisfying  the  damped  oscillator  equation 

p(t)  +  (iu^)^  +  y(t)  "  <*fi>2  0|^(u>,t)  »  (3.1) 

where  and  are,  respectively,  the  molecular  frequency  and  decay  rate  in 
the  absence  of  the  surface,  is  the  molecular  polarizability,  and  f(u),t)  = 
S(w)e  la)t  is  the  external  field  at  the  site  of  the  molecule.  Let  us  write 
this  field  in  the  form 

I(w)  =  iL  +  2  +  I  +1.  =  [1  +  A(w)  ] -In(u))  +  G(u))-p(w)  ,(3.2) 

u  r  sp  im  u 

where  each  term  stands  for  the  incident,  reflected,  surface  plasmon  and 
image  fields,  respectively,  and  the  coefficients  A(w)  and  G(oj)  are  in 
general  tensors.  A  Fourier  analysis  of  Eq.(3.1)  then  admits  the  solution 


m(oj)  = 
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where  n  is  the  unit  vector  of  the  direction  of  the  molecular  dipole  and  we 
M 

have  assumed  e  1U)t.  Note  that  and  are,  respectively,  the 

"surface-modified"  frequency  and  decay  rate  which  are  to  be  determined  by 
the  image  field  G(w) • p(w)  in  Eq.(3.2)  (Gersten  et  al.,  1985).  Using 
Eq.(3.3),  the  Poynting  flux  absorbed  by  the  molecular  system  has  been 
calculated  (Gersten  et  al.,  1985),  and  together  with  the  result  for  the 
incident  flux  this  leads  to  the  absorption  cross  section 


o(uj)  =  2irca2ar.  I  n  -[l+A]-n_|2 
0 1  u  O' 


(3.4) 


where  n^  =  Eq/eq,  a  is  the  fine  structure  constant  and  a^  the  Bohr  radius. 

It  is  clear  from  Eq.  (3.4)  that  o(oj)  exhibits  the  general  Lorentzian 
structure  and  that  surface  effects  enter  into  the  absorption/dissociation 
process  via  the  terms  A ( i u) ,  and  y^  .  Since  under  most  circumstances  the 

change  of  the  molecular  frequency  due  to  the  presence  of  the  image  field  is 
almost  completely  negligible  (Chance  et  al.,  1978),  we  shall  assume 
in  most  of  our  discussion  below.  It  then  becomes  clear  from  Eq.(3.4)  that 
two  competing  mechanisms  exist  for  such  processes  due  to  the  presence  of  the 
surface,  namely,  the  enhanced  local  field  (A(w)),  which  tends  to  increase 
the  rate  of  absorption/dissociation  processes,  and  the  increased  decay  rate 
(yu),  which  tends  to  suppress  such  processes  (Garoff  et  al.,  1982;  Gersten 
et  al.,  1985).  Because  of  this,  the  molecule-surface  separation  dependence 


(d)  for  this  kind  of  processes  is  very  different  from  that  for  ordinary  SERS 
(Murray,  1982).  For  the  latter  where  only  the  enhanced  SP  field  plays  an 
important  role,  the  enhanced  SERS  cross  section  is  expected  to  decrease 
monotonically  as  d  increases.  Because  of  this  difference  in  distance 
dependence,  we  shall  see  that  the  concepts  of  the  critical  and  optimal 
distance  mentioned  above  should  be  introduced  into  the  description  of  the 
photodissociation  process.  Before  we  present  the  results  from  model 
calculations,  however,  we  shall  summarize  some  work  dealing  with  the 

quantities  A(oj)  and  yw(oj). 

M 

The  problem  of  surface  electromagnetic  (EM)  fields  has  attracted  much 
attention  in  the  past  two  decades.  A  large  number  of  investigations  have 
been  made  in  order  to  determine  the  extent  to  which  the  classical 
macroscopic  Maxwell  theory  may  be  applicable.  This  covers  the  cases  for 
different  polarizations  (s/p)  of  the  incident  light,  different  kinds  of 
dielectric  (metallic/nonmetallic)  surfaces  and  different  surface 
morphologies  (flat/rough).  We  refer  to  the  recent  comprehensive  review 
article  on  this  subject  by  Feibelman  (1982),  particularly  his  detailed 
discussion  on  the  nonlocal  EM  theory  for  the  description  of  surface  optics. 
Here  we  shall  only  outline  briefly  the  case  for  a  rough  metallic  surface, 
with  special  attention  given  to  the  surface  plasmon  mode. 

Theoretical  modelling  of  surface  roughness  falls  into  the  following 
three  catagories: 

(i)  Periodic  roughness:  This  class  of  regular  roughness  is  normally 
represented  by  a  grating.  For  the  case  of  a  shallow  sinusoidal 
grating,  the  surface  fields  have  been  calculated  using  the 

Rayleigh  method  and  a  phenomenological  approach  (Marvin  et  al.,  1975; 
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Maradudin  et  al.,  1975;  Equiluz  et  al.,  1983;  Maradudin,  1982).  The 

case  for  a  deep  grating  (£,/A  >  0.072)  has  also  been  treated  using  the 

8 

extinction  theorem  (Toigo  et  al.,  1977)  and  the  dressed  Rayleigh 
expansion  approach  (Agassi  et  al.,  1986). 

(ii)  Random  roughness;  This  kind  of  surface  is  specified  by  a  spatial 
correlation  function  <C( r-r')>,  where  very  often  a  Gaussian 
distribution  of  roughness  is  assumed  so  that  <£(r-r’ )>  - 

exp(  |  r -r 1  |  ^/a*" )  with  both  coordinates  r  and  r'  lying  on  the  surface. 

The  surface  EM  fields  for  this  case  have  been  treated  using  a  Green 
function  approach,  mainly  by  the  Irvine  group  (Maradudin  et  al.,  1975; 
Equiluz  et  al.,  1983  Maradudin,  1982). 

(iii)  Island  surfaces:  This  kind  of  surface  is  usually  modelled  as  a 
collection  of  spheroids  which  can  be  spheres  or  ellipsoids, 
eitherregularly  (periodically)  or  randomly  distributed.  Normally  the 
case  for  one  spheroid  is  worked  out,  and  the  final  result  is  then 
obtained  by  the  principle  of  superposition  of  the  contributions  from 
the  other  spheroids.  The  surface  EM  fields  for  this  case  haveattracted 
great  interest  recently  (Meier  et  al.,  1985;  Cline  et  al.,  1986). 

As  for  our  application  here  to  the  dissociation  process  described  by 
Eq.(3.4),  we  shall  assume  the  simplest  geometry  to  illustrate  the  ideas  of 
the  various  competing  mechanisms  and  the  various  distance  concepts. 
Specifically,  we  shall  consider  the  rough  surface  as  a  shallow  sinusoidal 
grating  (see  Fig.  4).  As  mentioned  above,  the  surface  EM  fields  have  been 
worked  out  by  a  phenomenological  approach  with  the  application  of  the 
Rayleigh  hypothesis  (Marvin  et  al.,  1975;  Jha  et  al.,  1980).  Assuming  a 


notation  of  E  =  ")  with  Em  being  the  component  on  the  xy-plane  (i.e., 
z  *' 

along  the  surface)  and  the  z-component,  the  field-amplification  factor 
A( o.’ )  for  p-polarized  incident  laser  light  can  be  written  as  (Leung  et  al., 
1986) 


A(oj)  = 
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where  kg  =  ^  sin0  +  g,  g  =  2tt / with  being  the  spatial  period  of  the 

grating,  rg  =  (kg  -  =  ~  cos0  ,  k£  =  ^  sin0  ,  and  0  is  the  angle 

of  incidence.  The  quantities  R  and  S  in  Eq.(3.5)  ,  originating  respectively 
from  the  reflected  and  surface  plasmon  fields,  are  given  as 
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where  p 


2 


2 

~  c,  £  is  the  grating  amplitude  and  e  =  e(uj) 
c  8 


is  the  frequency-dependent  complex  dielectric  constant  of  the  metallic 
grating.  The  plasmon  resonance  condition  is  achieved  when 


(3.8) 


~  c.  (ui)  -  k  [e.  (w)  +  l  ]  =  0 
c2  1  8  1 

In  order  to  completely  investigate  the  surface  effects  entering  into 
Eq.(3.4),  we  must  know  the  induced  molecular  decay  rate  (y^)  at  the  surface. 
We  first  give  a  brief  review  on  this  subject  matter  for  the  case  of 
physisorbed  molecules. 

In  the  past  two  decades,  there  has  been  a  large  amount  of  effort  in  the 
study  of  the  lifetime  of  an  excited  molecule  in  the  vicinity  of  a 
surface, both  theoretically  (Kuhn,  1970;  Chance  et  al.,  1978)  and 
experimentally  (Drehage,  1974;  Rossetti  et  al.,  1980,  1982).  The 
theoretical  work  has  included  classical  reflected  field  theory  (Kuhn, 

1970)  and  the  energy  flux  method  (Chance  et  al.,  1978)  for  both  perfectly 
or  partially  reflecting  flat  surfaces  (Philpott,  1975).  The  classical 
approaches  have  also  been  extended  to  the  case  of  a  rough  surface  including 
both  the  randomly  (Arias  et  al.,  1982)  and  the  periodically  (Leung  et  al., 
1986,  1987b)  roughened  cases.  These  later  extentions  have  been  based  on 
theories  of  image  potentials  for  a  point  charge  located  near  such  surfaces 
as  established  by  the  Irvine  group  (Rahman  et  al.,  1980a, b).  Among  these 
investigations,  the  general  conclusion  has  been  reached  that  the  molecular 
lifetime  is  in  general  shortened  due  to  the  induced  decay  rate  by  the 
reflected  field  from  the  surface,  and  surface  roughness  further  enhances 
such  decay  rate.  Furthermore,  large  dependences  on  the  orientation  as  well 
as  the  distance  of  the  molecule  from  the  surface  have  been  observed  (Chance 


To  illustrate  here  how  the  effect  of  the  induced  decay  rate  (y  )  enters 

n 

into  Eq.  (3.4),  we  shall  discuss  y^  for  a  shallow  sinusoidal  grating.  By 
applying  the  general  theory  of  Rahman  and  Maradudin  (1980a),  Rahman  and 
Mills  (1980b)  have  obtained  the  image  potential  for  a  point  charge  e  near  a 
shallow  sinusoidal  grating  as 


<J>(z) 


,  ,  4K  (gz) 

|  £  g{(e-l)[gK  (gz)  +  4,(82))  +  - I - },  (3.9) 

8  (e+1 )  8  0  z  1  z 


where  Kq  and  are  the  modified  Bessel  functions.  For  a  perpendicular 
dipole  located  at  (0,0,d),  we  can  calculate  the  image  dipole  field  per 
dipole  moment  from  Eq.  (3.9)  to  be  (Leung  et  al.,  1986) 


GR(w) 


=  1  £-1 
(e+1 ) 


u  8  , ...  J  4  g{«Ke'1)(«  +  d}  +  dj(gKI  +  “d} 


3K 

+  VC+l)(gK0  +  ^i)) 

d 


(3.10) 


We  remark  that  this  result  is  first  order  in  £  g,  in  contrast  to  that  for  a 

g 


randomly  rough  surface  which  has  a  lowest  order  result  in  (£^g)  (Arias  et 


al.,  1982).  Following  similar  steps  as  in  Chance  et  al.  (1978)  and  Arias  et 


al.  (1982),  we  obtain  finally  the  y^  on  a  shallow  sinusoidal  grating  as 


0 


In 


(3.11) 


where  q  is  the  quantum  yield  of  the  emitting  state.  G  (uj)  is  the 
corresponding  function  as  in  Eq.  (3.10)  for  a  flat  surface  and  is  given  by 
the  Sommerfeld  antenna  theory  as  (Chance  et  al.,  1978) 


-2£.<i  3 

G  (u)  =  -k  |  i  du  e  , 

J0  *1 


(3.12) 


2  1/2 

where  <3  =  kd,  =  ( i ^  -  eH^)/(2.2  +  e^)  >  2.^  =  -  i ( 1  -  u  )  and 
e2  -  -  i(c  -  u2)I/2. 

With  A(u)  in  Eq.  (3.5)  and  in  Eq.  (3.11)  substituted  into  Eq.  (3.4), 
we  have  illustrated  with  some  numerical  calculations  the  surface  effects  on 
the  direct  dissociation  of  molecules  at  4500  A  on  a  silver  (Ag)  grating 

-3 

with  a  roughness  parameter  £  /A  =  8  x  10  (Leung  et  al.,  1986).  Figure  5 
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shows  the  distortion  of  the  original  Lorentzian  line  profile  for  different 
molecule-surface  separations.  We  observe  that  the  SP  resonance  (at  about 
2.895  eV  of  photon  frequency)  introduces  a  sharp  edge  so  that  the 
asymmetrically  distorted  profile  shows  similar  behavior  as  for  a  Fano 
profile  (Fano,  1961).  Furthermore,  the  double-peak  feature  for  the  case  of 
a  spheroid  substrate  (Weitz  et  al.,  1983)  is  not  observed  here  since  we  have 
adjusted  the  SP  resonance  frequency  to  lie  close  to  that  of  the  molecular 
resonance.  Such  an  adjustment  is  physically  possible  since  one  can  vary  the 
many  parameters  such  as  the  angle  of  incidence  (9)  and  the  grating  period 
(Ag)  so  that  Eq.  (3.8)  is  fulfilled  for  u  =  u^.  Figure  6  shows  the 
enhancement  factor  versus  molecule-surface  distance  for  various  incident 
laser  frequencies  .  The  results  show  clearly  the  existence  of  the  critical 
distance  (d  ,  below  which  o/o_  <  1)  and  the  optimal  distance  (d  ,  at  which 


o/Cq  is  maximum).  We  observe  also  that  under  the  SP  resonance  condition, 
one  can  still  have  a  large  enhancement  even  at  distances  far  from  the 
surface.  In  general,  the  SP  resonance  plays  a  more  significant  role  then 
the  molecular  resonance  in  such  processes.  With  the  existence  of  these 
distances,  the  experiments  mentioned  in  the  beginning  of  this  section  may  be 
understood,  at  least  in  a  qualitative  manner.  Furthermore,  the  concept  of 
the  critical  distance  may  lead  to  practical  applications.  As  an  example,  we 
suggest  that  if  one  could  coat  the  metallic  surface  by  means  of  the  "fatty 
acid  monolayer  assembly  technique"  (Kuhn,  1968)  so  that  all  the  molecules 
are  kept  at  a  distance  above  the  critical  distance,  one  could  then  guarantee 
that  surface-assisted  dissociation  is  maintained.  Some  future  directions  in 
this  regard  are  mentioned  in  Sec.  6. 

4.  Line  Shapes 

Related  to  absorption/dissociation  but  in  itself  a  much  broader  class 
of  phenomena  is  that  of  line  profiles  for  molecules  in  the  vicinity  of  a 
surface.  The  study  of  such  "surface  distortions"  of  molecular  line  profiles 
can  give  invaluable  information  concerning  the  details  of  the  molecular 
processes  as  well  as  the  mechanisms  by  which  the  molecule  and  surface 
interact  with  each  other.  As  already  mentioned  in  Sec.  3,  for  example,  the 
observation  of  the  "double-peak  feature"  or  the  "sharp-edge  window"  in  the 
absorption/dissociation  spectrum  signifies  the  effects  of  the  surface 
plasmon  on  the  molecular  processes.  In  this  section,  we  shall  give  a  more 
comprehensive  review  of  this  class  of  phenomena. 

While  most  of  the  theoretical  work  has  been  devoted  to  physisorbed 


adspicies,  there  has  also  been  some  work  considering  chemisorbed  molecules 
(Metiu  et  ai.,  1978).  In  the  work  of  Metiu  and  Palke  (1978),  the  infrared 


spectroscopy  for  the  vibration  of  an  atom  chemisorbed  on  a  solid  is  studied 
assuming  a  coupling  between  the  adatom  and  the  lattice  phonons  of  the 
substrate.  The  analysis  shows  that  such  coupling  is  not  sufficient  to 
explain  the  large  line-broadening  observed  experimentally  (Nakata,  1976), 
and  possible  explanation  is  attributed  instead  to  the  coupling  to  the 
electronic  degrees  of  freedom  of  the  metallic  substrate.  Indeed,  from  our 
previous  presentation  in  Sec.  3,  we  have  found  in  a  phenomenological 
approach  that  the  SP  coupling  leads  to  an  appreciable  broadening  in  many 
cases  (Nitzan  et  al.,  1981;  Leung  et  al.,  1986)  for  physisorbed  systems.  It 
is  therefore  worthwhile  to  extend  this  previous  work  to  the  case  of 
chemisorbed  systems,  first  phenomenologically  and  then  with  a  more  thorough 
microscopic  formalism.  We  are  at  present  persuing  this  problem  in  our 
laboratory.  For  physisorbed  adspieces,  effects  of  the  coupling  between  the 
system  and  the  substrate  via  excitations  of  the  electron-hole  (eh)  pairs, 
phonons  and  photons  have  been  condisered  for  vibrational  line  shapes  of 
admolecules  (Gadzuk  et  al.,  1984).  Distinction  between  vibrational 
dephasing  and  relaxation  has  been  emphasised  in  this  latter  work.  Such 
vibrational  dynamics  has  also  been  analyzed  later  in  a  fully  quantum 
statistical  formalism  where  the  anharmonicity  of  the  vibration  is  emphasized 
(Huang  et  al.,  1985).  In  particular,  for  the  system  of  OH  on  SiO^,  the 
dephasing  rate  via  phonons  is  found  to  be  considerably  faster  than  the 
energy  relaxation  rate  (Hutchinson  et  al.,  1986).  In  addition,  the 
vibration-rotation  spectrum  for  physisorbed  HC1  on  Ar(lll)  surface  has  also 
been  investigated  in  a  semiclassical  "trajectory  approach"  with  a  collapse 
of  the  R  and  P  band  structure  at  low  rotational  energies  being  observed 
(Adams,  1986).  Furthermore,  aside  from  those  phenomenological  approaches 


which  we  mentioned  earlier  (Nitzan  et  al.,  1981;  Leung  et  al.,  1986),  the 
effects  of  the  roughness  of  surface  on  the  fluoresence  and  absorption 
spectra  for  a  molecular  dipole  near  a  randomly  rough  metallic  surface  have 
also  been  formulated  in  a  fully  quantum  mechanical  fashion.  It  was  found 
that  surface-roughness  in  this  case  induces  distortions  in  the  flat-surface 
Lorentzian  toward  a  Gaussian  line  shape  as  the  roughness  increases  (Agassi, 
1986).  However,  it  was  criticized  earlier  that  these  vibrational  spectra  of 
coupled  adsorbed  molecules  may  not  be  Lorentzian  to  start  with, even  in  the 
flat  surface  case  (Langreth,  1985;  Sorbello,  1985).  In  a  very  clever 
analysis,  Langreth  (1985)  has  shown  that  the  line  shape  for  an  isolated 
vibrational  mode  of  adspecies  on  a  metallic  surface  is  necessarily 
asymmetric  in  the  presence  of  the  damping,  where  the  profile  is  more  of  the 
Fano  type  (Fano,  1961)  than  the  Lorentz  type.  This  is  a  good  example 
showing  that  the  study  of  line  profiles  can  lead  to  an  understanding  of  the 
coupling  mechanism  between  the  molecule  and  the  surface. 

Aside  from  vibrational  and  rotational  spectra,  electronic  and 
desorption  spectral  line  shapes  are  also  of  great  interest.  In  particular, 
electronic  excitation  by  a  UV-visible  laser  leading  to  desorption  of  an 
admolecule  has  been  recently  studied  by  Lin  et  al  (1987a).  The  result  shows 
that  the  laser-stimulated  desorption  (LSD)  spectra  (LSD  yield  vs  laser 
frequency)  exhibit  the  well  known  Beutler-Fano  (Fano,  1961)  behavior  due  to 
the  configuration  interaction  between  the  molecular  excited  electronic 
states  and  the  desorption  continuum.  These  investigators  have  also  studied 
the  case  of  IR  LSD  spectra  (Lin  et  al.,  1987b),  where  the  same  Beutler-Fano 
line  shape  is  found  to  exist  due  to  the  coupling  between  the  intramolecular 
optically  pumped  vibrational  modes  and  the  desorption  continuum.  The 


effects  of  inhomogeneous  broadening  on  these  IR  LSD  line  shapes  have  also 

been  investigated  (Gortel  et  al.,  1986;  Lin  et  al.,  1987b).  Because  of  the 

heterogeneity  due  to  different  adsorption  sites  and  imperfections  in  the 

underlying  substrate.strong  asymmetry  and  double-peak  features  in  the 
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absorption  spectrum  have  been  reported.  In  case  when  optical  vibrational 
modes  are  involved,  the  LSD  spectra  can  be  analyzed  as  a  superposition  of 
Beutler-Fano  bands  (Lin  et  al.,  1987b). 

In  spite  of  the  numerous  studies  of  surface  effects  on  the  absorption 
spectra  of  admolecules,  it  appears,  in  all  the  previously  mentioned 
references,  that  the  free  line  profile  (i.e.,  the  absorption  line  profile 
for  a  free  molecule  in  the  absence  of  the  surface)  has  almost  always  been 
assumed  to  be  symmetric,  often  of  a  Lorentzian  type,  so  that  asymmetric 
distortions  are  brought  about  by  the  presence  of  the  surface.  There  remains 
therefore  one  whole  class  of  problem  previously  uninvestigated,  namely  that 
when  the  free  line  profile  is  already  asymmetric  in  nature.  This  would 
include,  for  example,  processes  like  autoionization  and  predissociation  in 
molecular  systems.  In  a  recent  work  (Leung  et  al.,  1987),  we  have  reported 
some  preliminary  results  in  this  direction  based  on  an  extention  of  our 
previous  work  in  the  distortions  of  the  Lorentzian  profile  on  top  of  a 
silver  grating  (Leung  et  al.,  1986),  which  we  now  describe  below. 

In  analogue  to  the  "driven  damped  oscillator  model"  ( Eq .(3.1))  which 
describes  a  free  Lorentzian  profile  in  the  absence  of  surface  effects,  we 
have  adopted  the  mechanical  model  recently  proposed  by  Sorbello  (1985)  to 
describe  the  asymmetric  Fano  (1961)  effect.  The  model  consists  of  the 
coupling  of  a  normal  mode  (u)q)  to  a  lossy  broad-band  system  such  as  a 
viscous  bath,  with  both  the  oscillator  and  the  bath  being  described  by  one 


degree  of  freedom.  The  line-shape  function  obtained  in  this  model  can  be 
reexpressed  in  the  form'" 


(q+e0)2 

Iq(uj)  =  - - —  oQ(w)  ,  (4.1) 

1+eo 


where  q  is  the  asymmetric  profile  index  depending  on  the  ratio  of  the 
driving  force  on  the  oscillator  to  that  on  the  bath,  and  Og(w)  is  the 
background  intensity  due  to  direct  excitation  of  the  bath.  The  reduced 
energy  variable  Eg  in  Eq.  (4.1)  is  expressed  as 

2(uj-cd  -A(jj0) 
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where  Au>g  and  Yq  are  the  level  shift  and  decay  rate  of  the  system  , 
respectively.  In  general,  Awg  <<  uig,  and  hence 


E 

O 
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with  Au>  =  w  —  uig. 

In  an  analogous  way  as  for  the  "distorted  Lorentzian  profile" 

(Eq  .(3.4)),  we  have  obtained  the  surface  distorted  Fano  profile  in  the  form 
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*  We  have  changed  some  of  Sorbello's  (1985)  original  notation.  Here 
we  denote  every  quantity  which  refers  to  the  free  molecule  case  in  the  absence 
of  the  surface  by  the  subscript  "0". 


where  now 


Here  all  quantities  without  a  subscript  'O'  refer  to  those  at  the  surface, 

and  we  have  neglected  the  effect  of  the  surface  on  q  and  Wq  (Leung  et  al. , 

1987).  We  have  also  assumed  the  substrate  to  be  a  shallow  sinusoidal 

grating,  and  hence  the  quantities  A(oj)  and  y  are  given  as  before  in 

Eqs.(3.5)  and  (3.11)  for  a  molecular  dipole  located  at  (0,0, d)  and 

perpendicular  to  the  surface.  With  some  reasonable  forms  for  Oq(u>),  we  have 

examined  the  surface  distortions  for  the  cases  of  certain  autoionizat ion  and 

predissociation  absorption  spectra,  where  results  are  shown  in  Figs.  7  and 

8,  respectively.  From  these  preliminary  results  basing  on  such  model 

studies,  we  observe  the  familiar  double-peak  features  and  the  broadening  of 

the  oringinal  profile  window  at  the  steep  edge  near  the  low-frequency  end, 

while  at  the  high  frequency  end  of  the  distorted  profile,  the  surface 

plasmon  resonance  leads  to  a  new  window.  Furthermore,  we  observe  that  at 

such  molecule-surface  distances  (d  -  500  A)  under  the  model  calculation 

conditions,  a  surface  enhancement  effect  is  in  general  observed,  implying 

that  enhanced  molecular  photo-predissociation  may  also  be  possible,  provided 

that  the  molecule  is  not  located  too  close  to  the  surface  and  the  resonant 

plasmon  field  decays  very  slowly  in  the  direction  normal  to  the  surface.  It 

has  been  argued  previously  that  due  to  the  "slow  nature"  (on  a  time  scale  of 
-9 

-10  s)  of  these  unimolecular  processes  such  as  predissociation,  the 
surface  enhancement  effects  can  hardly  be  effective  due  to  the  large  induced 
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molecular  decay  rate  caused  by  the  presence  of  the  surface  (Chuang,  1982). 
However,  as  seen  from  our  preliminary  model  study,  such  an  argument  is  not 
conclusive  since  there  are  many  parameters  (e.g.,  degree  of  the  roughness, 
angle  of  incidence,  etc.)  that  one  can  adjust,  so  that  the  molecule  can 
still  experience  a  large  enhanced  field  at  a  distance  relatively  far  from 
the  surface  at  which  the  induced  decay  rate  is  of  minor  consequence.  Thus, 
it  would  be  interesting  to  conduct  a  more  detailed  study  of  such  possible 
enhanced  "slow  processes".  The  practical  realization  of  these  ideas  in 
actual  gas-surface  systems  might  lead  to  the  development  of  a  new  kind  of 
photochemistry,  namely,  laser-assisted  heterogeneous  catalysis  (Lin  et  al., 
1984)  via  surface  predissociation  of  physisorbed  molecule,  so  that  the 
problem  of  desorption  of  the  reaction  products  from  the  surface  can  almost 
be  completely  ignored.  We  should  also  add  that  possible  enhanced  molecular 
predissociation  has  been  examined  recently  by  means  of  yet  a  different 
mechanism,  namely,  surface  magnetic  field/laser  synergistic  effects 
(Bhattacharyya  et  al.,  1981,  1982).  In  contrast  to  the  SP-  enhanced  pumping 
rate  of  the  absorption/dissociation  process  discussed  here,  this  latter 
mechanism  creates  its  own  dissociation  channels  via  crossings  between  the 
ground  and  continuum  electronic  potential  surfaces  due  to  the  photon-dressed 
effect  associated  with  the  ground  state,  and  the  Zeeman  splitting  of  the 
continuum  levels. 

5.  Laser-Induced  Deposition 

A  relatively  new  field  has  emerged  in  which  a  low-power  laser  is  used 

to  enhance  deposition  of  a  metal  on  a  surface,  generally  semiconductor.  The 

2 

archetypal  experiment  involves  a  weak  257  nm  cw  laser  (1-10  W/cm*") 


irradiating  a  SiO^  surface  under  an  organometallic  gas  (CdlCH^^  at  1  torr 
pressure  (Brueck  et  al.,  1982),  for  example.  A  closely  related  subject,  but 
one  which  is  in  many  ways  similar,  is  the  topic  of  laser- induced  damage  at 
surfaces.  This  has  been  studied  much  more  extensively  (Fauchet  et  al. , 

1982;  Sipe  et  al.,  1983).  We  shall  find  that  the  deposition  problem  is 
simpler  in  that  the  relevant  chemistry  appears  to  take  place  well  above  the 
surface,  and  hence  the  Rayleigh  expansion  can  be  used  to  calculate  the  field 
strengths . 

The  common  thread  between  these  two  problems  is  the  interaction  between 
the  laser  and  the  surface.  Here  primarily  two  effects  predominate:  laser 
heating  and  the  excitation  of  surface  plasmons.  Sometimes  both  effects  are 
important,  though  in  our  deposition  example,  the  surface  heating  amounts  to 
only  about  1  K.  In  general,  however,  given  the  interaction  between  the 
surface  and  laser,  several  mechanisms  can  occur: 

1)  The  laser  can  melt  the  surface,  and  then  the  plasmon  or  diffusive 
motion  along  the  surface  creates  a  standing  wave  in  the  molten 
substance.  Upon  cooling  this  pattern  is  frozen  into  the  surface 
(Fauchet  et  al.,  1982;  Ehrlich  et  al.,  1983;  Osgood  et  al.,  1982;  Chen 
et  al. ,  1985) . 

2)  The  laser  can  heat  an  adsorbed  layer  (as  opposed  to  the  true  surface) 
and  cause  rearrangement.  This  has  been  observed  with  copper  deposition 
on  silica  surfaces  (Moylan  et  al.,  1986;  Ehrlich  et  al.,  1983;  Osgood 
et  al.,  1982;  Chen  et  al.,  1985). 

3)  The  laser  can  affect  the  adsorption  process  itself.  This  can  happen 
either  by  the  enhancement  of  the  adsorption  process  because  of  the 


laser  (we  are  unaware  of  any  experimental  evidence  for  this  occuring), 
or  the  laser  can  enhance  the  dissociation  of  the  molecular  precursors 
to  adsorption.  Brueck  and  Ehrlich's  experiment  appears  to  be  an 
example  of  this  latter  phenomenon. 

Sipe,  Young,  Preston  and  van  Driel  (1983)  have  written  extensively  on 
the  first  process.  To  account  for  this  phenomena,  it  is  essential  to  know 
the  field  intensities  at  the  surface.  As  mentioned  in  Sec.  1,  this  is  a 
very  complicated  problem.  They  have  calculated  the  energy  deposition  (due 
to  the  laser)  just  below  the  surface  and  then  used  a  variational  calculation 
to  compute  the  local  field  strengths  near  the  surface.  Moylan,  Baum  and 
Jones  (1986)  have  discussed  the  second  possibility  and  have  related  the 
grating  growth  rate  to  the  focal  point  size  of  the  laser,  along  with  the 
diffusion  coefficient  of  the  adsorbed  layer.  They  are  then  able  to  generate 
the  kinetic  equations  for  motion  along  the  surface. 

The  primary  subject  of  this  section  will  be  to  discuss  the  third 
process,  namely,  we  shall  try  to  account  for  the  observations  of  Brueck  and 
Ehrlich  (1982).  These  results  are  more  fully  described  in  Jelski  and  George 
(1987).  Put  briefly,  the  laser,  in  this  case  weak  so  as  to  minimize  surface 
heating,  causes  a  resonant  plasmon  due  to  coupling  with  the  surface 
roughness.  This  plasmon  enhances  the  dissociation  of  the  inorganic 
precursor  above  the  surface,  but  since  the  field  strength  is  spatially 
periodic,  this  leads  to  periodic  deposition. 

A  theoretical  treatment  of  this  process  requires  two  steps.  In  the 
first  step,  we  need  a  simple  method  to  calculate  the  field  strength;  the 
difficulties  involved  here  have  been  discussed  in  Sec.  2.  We  note  here  that 


the  problem  is  enormously  simplified  by  the  fact  that  dissociation  takes 
place  well  above  the  surface  (see  Sec.  3).  This  permits  use  of  the  Rayleigh 
expansion,  which  is  exact  above  the  selvedge  region.  Secondly,  we  want  a 
fast  calculation  of  the  field  strength  since  our  problem  requires  this 
strength  as  a  function  of  grating  wavenumber  and  grating  height.  The 
approximation  of  Jelski  and  George  (1987)  described  in  Sec.  2  seems  most 
appropriate  to  our  problem  since,  in  every  case,  we  require  the  intensity  of 
only  the  resonant  frequency.  A  general  reflectivity  calculation  is  not 
called  for. 

Thus,  stated  explicitly,  our  method  is  as  follows: 

1)  We  assume  (see  Jelski  et  al.  (1987b)  for  a  more  complete  discussion) 
that  the  enhanced  deposition  is  due  to  the  fact  that  the  plasmon  field 
enhances  the  dissociation  of  the  adspecies  precursor. 

2)  This  dissociation  takes  place  well  above  the  surface  (see  Sec.  3),  thus 
permitting  the  use  of  the  Rayleigh  expansion.  We  use  the  approximation 
of  Jelski  and  George  to  calculate  this. 

3)  To  simplify  the  calculation,  we  assume  that  the  dissociation  cross- 
section  of  the  precursor  is  proportional  to  the  field  intensity  and 
that  the  deposition  rate  is  proportional  to  the  dissociation  rate.  The 
first  approximation  breaks  down  for  strong  fields,  whereas  the  second 
fails  for  high  pressures. 

4)  Finally,  we  assume  that  the  final  grating  profile  is  sinusoidal.  This 
will  turn  out  to  be  approximately  accurate,  and  we  shall  see  from  the 
calculation  that  it  can  be  corrected,  at  the  cost  of  considerably  more 
complexity. 


Our  model  will  be  the  experiment  of  Brueck  and  Ehrlich  (1982).  We 

shall  assume  a  cadmium  thin  layer  with  a  dielectric  constant  of  -2.5  +  1.3i 

Dimethyl  cadmium  gas  at  ~  1  torr  pressure  and  at  room  temperature  is  above 

2 

the  surface.  A  257  nm  laser  at  10  W/cm  is  irradiating  the  system.  The 
dissociation  cross-section  of 


Cd(CH3)2  +  hv  -  Cd  +  2CH3 


(5.1) 


will  be  denoted  by  o^,  and  hence 


^dT1  '  »kEk  5  V(k)  • 


(5.2) 


where  E  is  the  field  strength  and  I  is  the  intensity.  The  subscript  k 
refers  to  the  wavenumber  of  the  incident  light. 

Before  going  further,  we  note  that  cadmium  is  a  lossy  substance  with  a 
large  imaginary  part  to  the  dielectric  constant.  Hence  the  flat  surface 
dispersion  relation,  given  by  Eqs.  (2.30)  or  (3.8),  is  inappropriate  since 
it  assumes  that  is  small.  We  thus  need  to  reconsider  the  flat  surface 
case  when  cannot  be  so  considered.  From  Eq.  (2.29),  and  noting  that 
0,  to  first  order  in  we  get 


ylEi  i(a0gl  "  k0kl)(eal  ‘  Pl)5lEi 


ek^  -  k^(l+e) 


(5.3) 


\.  '  V  %  ■ 


This  is  the  shallow  grating  limit  result  previously  derived  (Maradudin, 
1982;  Jha  et  al.,  1980),  which  is  identical,  except  for  notation,  to  Eq. 
(3.7).  Since  this  is  an  expansion  only  to  first  order  in  £^,  it  is  valid 
only  for  shallow  gratings.  Let  us  denote  the  denominator  by  R.  Then  if 


R  =  ek2  -  k2( 1+e)  =  0 


(5.4) 


Aj  becomes  large  and  resonant.  This  leads  to  the  resonance  conditi 


on 


.2  i2  1 

ki  = k  r+r  • 


(5.5) 


1 


where  e^  is  the  real  part  of  z.  Equation  (5.5)  is  essentially  the  condition 

expressed  in  Eq.  (3.8)  and  is  valid  if  is  small.  If  that  is  not  the 

2 

case,  then  Eq.  (5.4)  will  never  hold  and  instead  we  have  to  minimize  R  . 

This  yields 


| R | 2  =  [k2( 1+e. )  -  k2e  ]2  +  [c  (k2  -  k2)]2 
a  1  1  2  a 


(5.6) 


or  upon  setting  the  derivative  equal  to  zero. 


2  k2[e1(l+e1)  -  e2] 

k  - - 


8  (1+e, )2  +  e2 


(5.7) 


'V 


Equation  (5.7)  reduces  to  Eq.  (5.4)  as  -*  0. 


As  the  grating  becomes  deeper,  R  also  depends  on  powers  of  £  .  But  in 

m 

general,  we  can  express  the  resonance  condition  by 


dk. 


(5.8) 


For  very  deep  gratings  this  may  not  be  quite  sufficient,  as  we  would  also 
need  to  account  for  the  behavior  of  the  numerator,  although  for  our  purposes 
here  it  is  sufficient. 

Figure  9  shows  the  resonance  factor,  |R|,  as  a  function  of  grating 
wavenumber  for  the  model  cadmium  surface.  For  the  shallow  grating  (1  niru 
the  resonance  is  well  defined.  As  the  grating  deepens  it  becomes  apparent 
that  the  resonance  is  less  and  less  pronounced. 

Figure  10  is  a  graph  of  R'  vs.  grating  wavenumber.  Here  two  things  can 
be  seen:  one  is  that  as  the  grating  gets  deeper  then  the  resonance 
frequency  (where  the  curves  cross  zero)  shifts  toward  higher  frequencies. 
This  precludes  any  grating  from  growing  sinusoidally,  for  as  the  grating 
gets  larger  the  wavelength  gets  shorter.  Nevertheless,  we  can  approximate 
our  grating  as  sinusoidal  since  the  shift  is  relatively  small. 

Secondly,  we  note  that  the  24  nm  curve  never  crosses  zero.  Hence  there 
is  no  resonance  at  all.  This  means  that  grating  growth  will  stop.  This 
does  not  imply  that  deposition  will  stop,  but  only  that  deposition  will 
occur  evenly  across  the  surface,  i.e.,  our  thin  layer  will  get  thicker. 

Hence  peaks  and  valleys  will  grow  identically.  We  can  already  predict  the 
maximum  grating  height,  namely  48  nm  (2  *  the  amplitude).  This  is 


0 


approximately  half  of  what  Brueck  and  Ehrlich  (1982)  observe.  This  error 
may  be  due  to  our  approximate  (and  frequency- independent )  form  for  the 
dielectric  constant,  our  ignoring  thin  layer  effects,  or  our  approximate 
form  for  the  Rayleigh  expansion. 

Having  determined  the  maximum  grating  height,  we  now  need  to  determine 
the  rate  of  grating  growth.  Hence  we  must  ask  the  question  "How  much  more 
resonant  is  a  given  frequency  than  its  neighbors?"  In  a  word,  we  want  to 
know  the  concavity  of  the  curves  in  Fig.  9.  Let  us  define 


!  <X2>k  -k 
3  kg  6  r 


(5.9) 


where  k^  is  the  resonant  frequency  determined  by  Eq.  (5.8).  As  the 
resonance  disappears  the  concavity  goes  to  zero. 

Let  us  assume  a  sinusoidal  grating  of  amplitude  £,  with  N  as  the  number 
of  cadmium  atoms  that  adsorb  at  the  peak.  Then  we  can  write 


d£  -  dN 
dt  dN  dt 


(5.10) 


where 


Mk 

dN  2py 


(5.11) 


M  is  the  atomic  mass  of  cadmium,  p  is  the  density  of  metallic  cadmium  and  y 
is  the  dimension  of  the  surface  (parallel  to  the  grooves). 
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The  rate  dN/dt  is  proportional  to  right-hand  side  of  Eq.  (5.2),  i.e., 
to  the  number  of  Cd  atoms  that  adsorb  on  the  surface.  We  will  denote  the 
sticking  coefficient  by  the  letter  a.  Then 

rlN 

g  -  aokI(k,r)G(kr,e)  (5.12) 

where  r  is  the  spatial  variation  of  the  intensity.  The  quantity  G(kr>£) 
weights  the  function  according  to  the  frequency  variation  of  the  resonance. 
Equation  (5.12)  can  be  integrated  to  yield  £(t).  Choosing  the  constants  to 
reflect  the  time  span  of  Brueck  and  Ehrlich's  experiment,  we  generate  the 
function  shown  in  Fig.  11. 

Thus  we  have  developed  a  simple  model  to  account  for  the  chemical  vapor 
deposition  phenomena.  We  have  assumed  that  the  ultimate  cause  of  this 
effect  is  the  laser-induced  plasmon  which  enhances  dissociation  above  the 
surface.  We  have  succeeded  in  qualitatively  reproducing  Brueck  and 
Ehrlich's  experiments.  A  more  sophisticated  approach  will  have  to  allow  for 
the  thin  layer  effect,  provide  a  more  accurate  function  for  the  dielectric 
constant,  relieve  the  condition  that  the  grating  be  sinusoidal,  and  perhaps 
use  the  full  Rayleigh  expansion  rather  than  the  approximate  form. 

6.  Conclusion 

In  this  paper,  we  have  presented  a  review  of  different  types  of 
photochemical  processes  at  rough  (metallic)  surfaces,  giving  attention  to 
the  work  done  by  others  as  well  as  by  our  group.  In  contrast  to  a  previous 
review  article  by  Harris  and  coworkers  (Goncher  et  al.,  1984)  which 
emphasizes  experimental  studies  in  photodissociation  processes,  our  present 


article  emphasizes  the  theoretical  aspects  and  incorporates  also  phenomena 
like  molecular  lineshapes  and  dposition  on  rough  surfaces.  Since  we  have 
been  following  the  classical  phenomenological  treatment  (Maxwell's  theory) 
for  the  surface  fields,  we  have  also  provided  the  mathematical  basis  (Sec. 

1)  that  clarifies  the  foundation  of  the  Rayleigh  expansion  which  has  been 
found  invaluable  in  the  treatment  of  surface  EM  fields. 

The  main  physical  mechanism  involved  in  all  the  processes  we  have 
discussed  has  been  the  surface  plasmon  effect.  This  has  turned  out  to  be 
the  case  since  we  have  assumed  physisorbed  systems  throughout.  For 
chemisorbed  systems,  bond- format  ion  effects  will  arise  between  the 
adnolecule  and  the  surface,  and  charge-transfer  processes  should  be  given 
equal  importance  as  the  SP  effect.  We  hope  to  generalize  our  previous  work 
in  this  direction  in  the  future.  Such  generalization  may  also  find 
applications  in  the  theoretical  exploration  of  possible  laser-assisted 
heterogenous  catalysis  as  well  as  deeper  understanding  of  the  deposition 
process  discussed  in  Sec.  4,  which  has  hitherto  been  treated  only 
phenomenologically.  Furthermore,  we  have  in  this  review  assumed  an  isolated 
molecule-surface  system.  An  actual  physical  situation  will  mor  generally 
involve  an  ensemble  of  molecules  in  the  vicinity  of  the  surface.  We 
therefore  also  plan  to  investigate  in  the  future  the  effects  of  the 
neighboring  admolecules  on  the  phenomena  reviewed  in  Secs.  3  to  5. 

Finally,  it  is  clear  in  this  review  on  surface  photochemistry  that  all 


the  analysis  is  based  on  the  linear  Maxwell's  theory.  Since  nonlinear 
optical  surface  processes  have  attracted  much  attention  recently  (Chen  et 
al.,  1981),  we  feel  that  the  time  is  ripe  to  look  into  other  possible 


photochemical  processes  exploiting  these  nonlinear  optical  techniques.  We 
have  begun  our  first  attempt  to  reinvestigate  the  above  processes  on  a 
phase-conjugated  surface,  replacing  the  ordinary  substrate  metallic  surface 
by  an  optically  nonlinear  substrate  such  as  GaAs(llO).  Some  interesting 
preliminary  results  have  been  obtained  with  respect  to  the  linewidth  and 
lifetimes  for  molecules  near  such  conjugated  surfaces  (Lin  et  al.,  1987). 

We  believe  that  it  would  be  very  worthwhile  to  invest  more  efforts  in  this 
direction  in  the  future. 
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run) 


R  (exact) 


R  (Weber's  method) 


Error  ( 7. ) 


0.855 

0.856 

0.117 

0.566 

0.570 

0.707 

0.163 

0.155 

4.91 

0.184 

0.158 

14.1 

0.298 

0.273 

8.39 

0.380 

0.363 

4.47 

0.422 

0.413 

2.13 

Percentage  error  of  Weber's  mode-coupling  perturbation 
calculation:  reflectivity  versus  height  over  the  range  of  grating 

heights  for  which  the  exact  theory  (undressed)  converges. 

Frequency  and  angle  (2.215  eV  and  32°)  were  chosen  so  that  that 
the  rainigap  region  of  the  surface  polariton  dispersion  curve  is 
being  probed.  The  surface  is  a  silver  sawtooth  grating  with  a 
wavelength  of  8000  A.  (This  table  is  reproduced  from  Weber  (1986) 
with  permission  of  the  author.) 


FIGURE  CAPTION’S 


Reflectivity  as  a  function  of  photon  energy  of  the  sawtooth  profile 
grating  for  two  values  of  the  grating  height  h.  The  solid  line  is  the 
exact  result  and  the  dashed  line  is  the  perturbation  theory  result.  In 
the  first  graph,  the  two  curves  are  coincident  to  within  graphical 
accuracy.  The  second  graph  is  done  for  the  maximum  grating  height  for 
with  the  exact  theory  converges.  (Reproduced  from  Weber  (1986)  with 
permission  of  the  author.) 

Plot  of  the  plasmon  field  strength  (A^)  calculated  from  Eq.  (2.26) 
using  the  dressed  Rayleigh  expansion  described  in  the  text.  This  is 
compared  with  the  exact  calculation.  The  dielectric  constant  used  is 
-2.5  +  1.3i,  the  incident  light  has  a  wavelength  of  257  nm,  and  the 
grating  wavenumber  is  taken  as  2.95  *  lO'7  m  *. 

Geometrical  arrangement  of  the  photodissociation  process. 

(a)  Sinusoidal  grating.  The  cross-hatched  area  represents  the  metal. 

(b)  Square-well  grating  showing  a  separation  into  three  layers,  one  of 
which  is  periodic  in  the  x-direction  and  two  of  which  are  uniform.  (c) 
Generalization  of  the  square-well  grating  in  which  there  are  three 
periodic  layers. 

Distortion  of  the  Lorentzian  line  profile  for  various  molecule-surface 
distances  for  a  perpendicular  molecular  dipole.  The  system  consists  of 
an  I2  molecule  on  a  silver  sinusoidal  grating.  We  refer  to  the  text 
for  numerical  data. 

Enhancement  factor  vs.  molecule-surface  distance  for  various  incident 

laser  frequencies,  (a)  E^  =  2.895  eV  (at  plasmon  resonance);  (b)  E^  = 

2.5  eV;  and  (c)  E  =  2.755  eV  (at  molecular  resonance).  Other 
Y 

parameters  are  the  same  as  in  Fig.  4. 


Distortion  of  the  Fano  profile  with  Og C )  a  constant  which  simulates 

certain  autoionization  processes.  The  profile  constants  are  ^  =  5  * 

10“,  Wq  =  1.6  *  104  and  q  =  -2.65.  Note  that  the  scales  for  I g C co )  and 

I(io)  are  different.  The  y-axis  quantities  are  in  arbitrary  units. 

2 

Distortion  of  the  Fano  profile  with  Oq(uj)  =  1/w  which  simulates 

certain  predissociation  processes.  The  profile  constants  are  =  4  * 

3  4 

10  ,  u»q  =  1.17  *  10  and  q  =  3.0.  The  y-axis  quantities  are  in 
arbitrary  units. 

Grating  wavenumber  vs.  the  value  of  the  resonance  factor  for  different 
grating  amplitudes. 

Value  of  R1  vs.  grating  wavenumber  for  different  grating  amplitudes. 
Note  the  frequency  shift  in  the  resonance  frequency  and  that  for  £  >  24 
nm  there  is  no  resonance  frequency. 

The  dynamics  of  grating  growth:  time  vs.  grating  amplitude. 
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